High frequency dynamics in liquid nickel: an IXS study 
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Owing to their large relatively thermal conductivity, peculiar, non-hydrodynamic features are 
expected to characterize the acoustic-like excitations observed in liquid metals. We report here an 
experimental study of collective modes in molten nickel, a case of exceptional geophysical interest for 
its relevance in Earth interior science. Our result shed light on previously reported contrasting evi- 
dences: in the explored energy-momentum region no deviation from the generalized hydrodynamic 
picture describing non conductive fluids are observed. Implications for high frequency transport 
properties in metallic fluids are discussed. 
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I. INTRODUCTION 

Half a century of inelastic neutron scattering experi- 
ments, recently complemented by similar investigations 
with X-rays, have clearly shown that acoustic like exci- 
tations can be sustained by simple liquids down to wave- 
lengths comparable to the mean interparticle distances 
and frequencies extending up to the THz region, In 
ordinary, non conductive fluids, such collective modes 
can be roughly described in terms of adiabatic sound 
waves of energy E = Ml propagating with sound ve- 
locity Cs = limQ_>o ^^Q^ ■ This is a well known, imme- 
diate consequence of the hydrodynamic treatment when 
the condition >> DtQ^ holds, where Dt is the ther- 
mal diffusion coefficient for the system under study and 
Q the exchanged momentum between the probe and the 
sample»ii^ A deeper analysis reveals that such acoustic 
waves are actually subject to relaxation processes related 
to the frequency dependent viscosity. The origin of such 
processes is twofold: the structural relaxation, respon- 
sible for dynamical arrest in those systems capable of 
supercooling, and a microscopic process which induces 
an additional damping of the sound waves due to the 
non-plane wave nature of the instantaneous vibrational 
eigenmodes. Both these processes manifest themselves 
at wavevectors Q w 1/10Qm, being Qm the principal 
maximum of the static structure factor. This region, un- 
fortunately, represents the technical limit for both neu- 
tron and X-ray inelastic spectroscopy, hindering the com- 
plete understanding of the underlying dynamical pro- 
cesses which are still controversially debated 

In liquid metals, an additional complication arises. 
The hydrodynamic condition ^{Q) >> DtQ^ breaks 
down in the range 0.1 ^ Q < 1 nm"-'^, an esti- 
mate obtained neglecting the Q dependence of the ther- 
mal conductivity. Consequently, an isothermal regime 
may be expected to occur above the fl{Q) ~ DtQ^ 
crossover, with the adiabatic limit attained only below 
this valuei^i^!^ Among the implications of such an isother- 



mal regime on the sound waves propagation, noteworthy 
features would be: 

• a reduced value of the sound velocity Ct = Cs/^/^, 
with 7 = CpjCvi the ratio of constant pressure to 
constant volume specific heats. 

• a different expression for the thermal contribution 



to the acoustic attenuation: 



instead of (7 — 



To address this issue, we present here an investigation of 
the high frequency dynamics in liquid nickel performed 
by Inelastic X-ray Scattering (IXS), a technique which 
has proved, since the early 90's, to be an invaluable tool 
for deepening our comprehension of the dynamics of sim- 
ple liquids on the microscopic scale, since it allows on 
one hand to overcome the kinematic limitations due to 
the lower incident energies of the neutrons and their spe- 
cific energy-momentum relation, and on the other hand 
provides direct access to the coherent cross section of the 
scattering process;^ 

In this respect, previous investigations performed by 
means of Inelastic Neutron Scattering (INS) and molecu- 
lar dynamics (MD) have shown contrasting resultsi^'^'^ 
In the former study the sound velocity attains indeed 
the isothermal value at the longest accessible wavelength 
{Q sa 8nm^^), while in the latter the numerically es- 
timated sound velocity is always larger than the adia- 
batic value. Recently an interesting study involving both 
quasi-elastic Neutron Scattering and MD has pointed out 
how the simulation data for the speed of sound of liquid 
Ni attain the adiabatic value cs at the lower accessed Q's 
(g«2nm-i)^ 

Liquid Ni, indeed, is characterized by the largest spe- 
cific heat ratio among monatomic liquids, a property 
which should emphasize any non hydrodynamic results. 
In addition, Ni is a system of paramount relevance in 
geophysical science due to its presence in the Earth's 
interiorjii Indeed, the fact that the Earth outer core is 
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mostly iron was established beyond reasonable doubt al- 
ready in the early 60 's, when it was confirmed that the 
density of the core was about 10% lower than the density 
of iron, and that the seismic parameter $ — K/ being 
K the bulk modulus and p the density, was higher than 
that of iron.^^ About 4% Ni is thought to be present in 
the core and although it does not appreciably change the 
density of liquid Fe, its presence should not be forgotten 
as phase diagrams of systems Fe-Ni-light elements may 
be significantly different from those of systems without 
■jsjj j^ji4 Among the most favorite candidates as lighter al- 
loying elements with Fe and Ni, sulfur is perhaps the most 
addressed oneJ^ In this respect, recently a high anoma- 
lous behavior of the ultrasonic sound velocity cs(T) and 
attenuation r(r) as functions of temperature was re- 
ported in the mixture 85%Fe-5%Ni-10%S, for temper- 
atures above melting (Tm =1650 K) up to 2000 K, at 
ambient pressure conditions. In fact, and contrary to 
the data of pure liquid metal components, the acoustic 
velocity is found to increase with temperature, as well as 
attenuation. A complete understanding of the underlying 
behavior of the alloy requires, in our opinion, to ascertain 
the nature of acoustic excitations in pure liquid Ni. 

II. THEORETICAL BACKGROUND: 
EXPECTED HYDRODYNAMIC BEHAVIOR FOR 
LIQUID METALS 

In an IXS experiment the double differential cross- 
section, which depends on the exchanged momentum Q 
and energy is proportional to the so called dynamic 
structure factor S{Q,u}), which in turn is the Fourier 
transform of the time dependent intermediate scattering 
function 

F{Q,t) ^l-Y^(e-'Q-meQ-^^''^y (1) 

Here N is the total number of particles constituting the 
system, and ri{t) the position of particle i at time t? 
In particular, the zero time value of the intermediate 
scattering function is directly related to the structural 
features of the system, being F{Q,Q) = S{Q), i.e. the 
static structure factor. In the hydrodynamic limit, for 
low values of Q, the dynamic structure factor displays 
three distinct peaks. A quasi elastic one, located at zero 
energy exchange, whose width Fe = DtQ^ is related to 
the thermal diffusion coefficient Dt, and two inelastic 
peaks - the so called Brillouin doublet - located at fre- 
quencies uj = zLcsQr^ with eg being the adiabatic speed 
of sound, and whose width depends mostly on kinematic 
viscosity, which is the main mechanism driving sound 
damping. On the basis of linear hydrodynamics, the dis- 
persive behavior, i.e. the dependence of the frequency uj 
of propagating collective modes on Q, is well known to 
display a transition between a linear adiabatic regime to 
a linear isothermal one, characterized by an isothermal 
speed of sound ct, at Q values such that ujTth ~ 1, where 



Tth is the characteristic decay time of thermal fluctua- 
tions {rth ~ In other words, at sufficiently 
high Q values, thermal fluctuations are expected to decay 
on a timescale much longer than the timescale of sound 
propagation, which now takes place in a thcrmalizcd en- 
vironment. If we look at the lineshape of the dynamic 
structure factor, the linewidth of the quasi-elastic line 
will increase, ultimately overwhelming the Brillouin dou- 
blet and causing a shift of the inelastic peaks position. 
The adiabatic and isothermal speed of sound, as we al- 
ready pointed out in section |T1 are related by the specific 
heat ratio 7, and we may call Qo ^ cs/Dt the value at 
which the adiabatic to isothermal transition takes place. 
In this scenario liquid metals constitute a particularly 
interesting class of simple liquids, because of their high 
thermal diffusivity with values of about ten times larger 
than cs, shifting the Q-range of the isothermal region 
down to values between 0.1 nm^^ and 3nm^^ provided 
that: 

1. the Q dependence of transport coefhcients can be 
neglected below at least 3 nm^^, as is the case en- 
countered in most liquid metals. 

2. the system interaction with radiation can be de- 
scribed within an effective single-component model, 
in which only the core electrons from the metal- 
lic ions interact with the electromagnetic radiation. 
However transport coefficients are supposed to take 
an effective value also carrying the net effect due to 
valence electrons. 

The second assumption has been on the basis of a recent 
debate focused on the interpretation of IXS data on liquid 
alkali metals and aluminum i^i^i^^ The main difficulty 
in ascertaining the existence of sound propagation with 
lower than adiabatic sound velocity value lies in the fact 
that 7 is for most metallic systems very close to unity, 
as can be seen from the specific heats ratios reported for 
a selection of liquid metals near their melting point in 
table [J For alkali metals 7 ~ 1 as well as for the great 
majority of liquid metals. Among the systems exhibiting 
a comparably high value of the specific heat ratio (Fe, Co 
and Ni), liquid nickel has been chosen, also in view of its 
lower - despite still very high - melting temperature. 

III. THE EXPERIMENT 

The experiment reported in this work was carried 
out at the high resolution Beamline ID16 of the Euro- 
pean Synchrotron Radiation Facility (Grenoble, Fr). The 
backscattering monochromator and analyzer crystals, op- 
erating at the (11,11,11) silicon reflections gave a total 
energy resolution of 1.5 meV, while energy scans were 
performed by varying the temperature of the monochro- 
mator with respect to that of the analyzer crystals. The 
sample, in the shape of a 100 pm thick foil, was 99.993% 
purity nickel, purchased by Rare Metallic Co. Ltd. 
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TABLE L Thermal properties of selected liquid metals near 
the melting point. 7, the ratio of specific heats, and Dt, the 
thermal diffusion coefficient, are reported at those temper- 
atures T, near the melting point, where data are available. 
For a more extensive summary of data on thermal as well as 
dynamical properties of liquid metals, see Ref. 3. 



Sample 


r[K] 


7 


Dt [nm^/ps] 


Li 


453 


16,17 18 

1.08— '—1.063— 


19 19 

19.1— ,20.3— 


Na 


371 


20 18 

1.12—1.091— 


19 

68.8— 


Mg 


923 


16 

1.29 


19 

39.8— 


Al 


933 


16 
1.4— 


19 

35.2— 


K 


336.7 


20 18 

1.11—1.102— 


19 

81.4— 


Fe 


1808 


16 

1.8— 


7.3ia 


Co 


1765 


16 

1.8- 




Ni 


1728 


16 

1.98— 


9.6iS 




1763 


6 

1.88- 


19 

42.1— 


Cu 


1356 


16 

1.33— 


Zn 


693 


20 16 

1.23—1.26— 


19 

15.7— 


Ga 


303 


1.08^^ 


11.6i^ 


Ge 


1253 


16.21 

1.18—'— 


19 

61.3— 


Rb 


312 


20 18 

1.13— ,1.09'^ 


Ag 


1233 


16 

1.32— 


19 

66.3— 


Sn 


505 


20 

1.11— 


19 

17.3— 


Cs 


302 




19 

44.6— 




308 


22 

1.102—, 


19 

40.4— 


An 


1336 


16 

1.28— 


Hg 


293 


23 

1.14— 


19 

4.41— 


Pb 


623 


16.17 

1.19—'— 


19 

9.89— 
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FIG. 1: (Color online) a: sound dispersion for liquid Ni: 
IXS (□) from the present work; an INS investigation by 
Bermejo et. al^ (A); a MD simulation from Ruiz-Martm et 
al. ( }^. Adiabatic ( — ) and isothermal ( ) sound dis- 
persions are also reported, b: Speed of sound reported from 
the present work (□) and from INS- (a) and MD (- • 
and oi) investigations. The adiabatic, cs, and isothermal, 
ct, speed of sound are displayed (• • ■). 



IV. RESULTS AND DISCUSSION 



(Japan), and was hold in a sapphire ceU obtained from a 
single monocrystali^ The use of sapphire prevented the 
sample interaction with container, still providing a good 
cell transmission, with a total cell thickness of 500 /im 
seen by the scattered beam. The cell was lodged in a 
molybdenum holder, and heated up by means of a prop- 
erly isolated tungsten resistance. The nickel absorption 
length at an energy of the incoming beam of 21 KeV, is 
about 50 ^m, thus in our condition a 15% sample trans- 
mission was expected. A selection of IXS spectra from 
liquid Ni at 1767 K is reported in Fig. [T] as open circles 
at the lowest accessed Q values. Dotted lines report the 
experimental resolution. Phonon modes from the sap- 
phire cell (speed of sound of about llkms"^) are well 
recognizable up to 3nm~^, however they are well sepa- 
rated from inelastic features from the sample at Q greater 
than 1 nm"""^. 



Experimental data have been analyzed with a damped 
harmonic oscillator (DHO) function centered at fre- 
quency ujl^ and of width modelling inelastic con- 
tributions from the metal. The elastic features have 
been represented by a Lorentzian of full width at half 
maximum (FWHM) Tc- Such an approximation is ex- 
pected to work well for not too high Q's, i.e. until there 
is a clear separation between elastic and inelastic fea- 
tures, as is the present case for IXS measurements. At 
Q values approaching Qm, corresponding to the static 
structure factor maximum, one has to resort to a model 
based on a so-called extended heat mode and two ex- 
tended sound modes, also incorporating the frequency 
sum rules.— The elastic contribution from the cell, sup- 
posed to be much narrower than the experimental resolu- 
tion, has been modelled with a Dirac delta function cen- 
tered at zero frequency. Thus the model for the S{Q,lo) 
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Q [nm"^] 

FIG. 2: (Color online) a: sound dispersion for liquid Ni: 
IXS (□) from the present work; an INS investigation by 
Bermejo et. alfi (A); a MD simulation from Ruiz-Martin et 
al. ( }^. Adiabatic ( — ) and isothermal ( ) sound dis- 
persions are also reported, b: Speed of sound reported from 
the present work (□) and from INS*' (a) and MD (- • 
and c^) investigations. The adiabatic, cs, and isothermal, 
Ct, speed of sound are displayed (■ • ■). 
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FIG. 3: (Color online) a: sound damping from the present 

work (□) and from INS** (a). The dashed line ( ) is 

a quadratic best fit to the IXS data, b: the quasielastic 
linewidth obtained in the present investigation of liquid Ni 
(□) and from a recent INS investigation (from the coherent 
dynamic structure factor (♦). The expected hydrodynamic 

behavior is also reported ( ), from the knowledge of the 

thermal diffusion coefficient 



can be written as 



Iq{Lo) ^ IcAQ)5{oj) + IcAQ) 



rc(Q) 



-Il{Q) 



(2) 



Additionally, in order to properly reproduce the experi- 
mental spectra, the above expression has to be modified 
such as to be compliant with the detailed balance condi- 
tion: 



(3) 



Before comparison with IXS data, the model is convo- 
luted with the experimental resolution. This procedure 
results in the curves displayed as full lines in Fig. [TJ Fig. 
[2] reports the dispersion (a) and sound velocity (b) ob- 
tained for liquid Ni from the DHO model parameter oj^ 
at different Q values (□), while Fig. [31a reports sound 



damping as derived from the DHO parameter (□). 
Sound speed and damping data are also reported from 
the cited INS experiment (a)^, as well as from MD ex- 
periments (ol, — • — ^). On one hand our results (□) for 
the acoustic excitations frequency dependence on Q are 
in very good agreement with the data obtained by Ruiz- 
Martin et ali^ ( ) in the most recent MD simulation 

on the subject, while on the other hand the experimen- 
tal INS data (A) are qualitatively different from both 
sets of data along the momentum range investigated.— 
This discrepancy is particularly evident in the region of 
Q <12 nm~^, where INS estimated sound speed reaches 
the isothermal value ct (Fig. [2]). This could be ascribed 
to the kinematic limitations to the energy accessible win- 
dow at the lowest Q's which prevented the observation of 
the S{Q,uj) tail, and thus a reliable estimate of the DHO 
model parametersi^ 

In order to address the last issue, in Fig. [4] we report a 
comparison of IXS (o) and INS (A,^) dynamic structure 
factors a.t Q — 8nm~^, corresponding to the lowest 
accessed Q-point for neutrons. The inelastic contribu- 
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tion to the IXS spectrum is emphasized by reporting 
the DHO function derived from the model (Eqs. 
and dH])). Beside the quasielastic peak, reflecting also 
the effect of incoherent scattering (which is not present 
in IXS in the case of monatomic systems), INS data 
show on the anti-Stokes side the presence of an inelastic 
peak around -16 meV. However, although the detailed 
balance condition would imply an even more pronounced 
peak on the Stokes side, its symmetric counterpart 
cannot be observed. Last but not least, the energy of 
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FIG. 4: Comparison of IXS (o) and INS (Ap dynamic struc- 
ture factors for liquid Ni at Q = 8nm~^. Data have been 
normalized to a constant factor in order to plot them on the 
same scale. The DHO function used to model IXS inelastic 
features is also reported ( — , see also Fig. [TJ. The exchanged 
energy range covered by INS is narrower than that accessed 
by IXS, due to the kinematics of the scattering process for 
neutrons. 

propagating excitations found by means of IXS (22.6 
meV) lies just at the edge of the accessed energy range 
for INS (-20 to 20 meV; see Fig. |4]), corresponding to 
the energy and scattering angle configurations used in 
this experiment. Summing up, while at sufficiently high 
Q's the sound dispersion of liquid Ni lies well inside 



the region of the momentum-energy range covered by 
INS, at Q values approaching wavevectors as low as 
8nm^^ the reliability of the dispersion curve extracted 
by INS is questionable. In Fig. [3]b, we report data 
for the elastic linewidth, compared with the expected 
hydrodynamic behavior, in which case the width of the 
quasielastic line is entirely due to thermal relaxation 

Tc = DtQ^ { ). Experimental INS values for Fc 

are also shown from a recent work (♦).^*' The present 
set of data (Fig. [31b) seems to unsupport the hypothesis 
of a weak Q dependence of transport coefficients 4 
Under this circumstance, the crossover to an isother- 
mal regime in liquid Ni may occur at values well below 
1-^2 nm^^, thus not accessible in the present experiment. 

V. CONCLUSIONS 



In conclusion, we performed an IXS experiment on liq- 
uid Ni at 1767 K, showing the ability of such system 
to sustain sound propagation over wavelengths compa- 
rable to the typical interparticle distance. Despite the 
indication from a previous INS investigation^ suggest- 
ing the tendency of sound speed to attain the isother- 
mal value at Q lower than 8nm~^, we found the evi- 
dence of an adiabatic dynamical regime holding at Q as 
low as 2nm~^, thus confirming MD results for the mi- 
croscopic dynamics of this systemiiiiS The discrepancy, 
observed for 8 < Q < 12 nm~^ between INS data on 
one hand and results from IXS and MD on the other 
hand, is probably related to the limitations imposed to 
the energy-momentum accessible region by the kinemat- 
ics of the scattering process for neutrons. Even though an 
adiabatic regime still holds at the lowest accessed Q for 
liquid Ni, the weak Q dependence of the quasielastic line- 
width clearly suggests that a generalized hydrodynamic 
picture should be invoked. The present result calls for 
further investigations with higher signal to noise ratio 
and at different temperatures to unambiguously ascer- 
tain the existence of an intermediate isothermal regime 
in liquid Ni in the explored momentum region. 
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